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Abstract 

The objective of this work is to reshape a nacelle to achieve a specified nacelle pressure distribution. 
The nacelle may be either isolated or installed on an airplane. There are no restrictions on the attitude 
(toe, incidence, and roll) and position of the nacelle. The design algorithm is coupled to two different 
multi-block 3-D Navier Stokes flow solvers. The coupling between design and analysis is automated to 
the point where the design proceeds with minimal user input. 

Nomenclature 

C - Local nacelle or wing chord. 

CDISC - Constrained DISC. 

Cp - Coefficient of pressure. 

DISC - Direct Iteration Surface Curvature. 

0(..) - Order of magnitude. 

R - Local nacelle radius. 

X - Streamwise coordinate. 

theta - Angle running from the nacelle crown through the keel (outboard side) and back up to the crown 
(inboard side). 

Introduction 

Given the harsh reality of placing capital intensive products in the modern marketplace, it is an ongoing 
quest to improve the operational efficiency and cost effectiveness of any commercial airplane. The 
emergence of the current generation of large turbofans, in the 80,000 to 100,000 pound thrust class, 
presents anew the opportunity to optimize or reduce the drag associated with the outer surface of the 
engine nacelle. These large nacelles could account for a good portion of airplane drag. 

The problem of external cowl design may be broadly placed in two separate, yet complementary 
categories: (1) the design of an isolated cowl immersed in a uniform flow; and, (2) the design of an 
installed cowl with the attendant non-uniform flow. The former should enable the designer to derive a 
cowl geometry starting from a fairly rudimentary first guess. The resulting design will be inherently 
symmetrical about the vertical plane, thus allowing for a universal nacelle that can be installed on either 
side of the airplane. On the other hand, the inverse design of an installed nacelle will always result in a 
nacelle that is non- symmetric about the vertical plane, unless the designer can invoke special options to 
force symmetry. 



Thus, a tool for the inverse design of nacelles must be flexible enough for a designer to be able to 
produce an initial design, then study changes required to the proposed design because of the presence of 
the airplane flowfield, and finally, perform trade studies of the pros and cons of nacelle symmetry. For 
example, the cruise benefit gained from achieving laminar flow over a non-symmetric nacelle might be 
significant enough to outweigh other issues. 

Method 

The design method used in the codes is the CDISC inverse design algorithml. CDISC is an extension of 
the basic DISC method of Campbell2&3. At the very crux, DISC is a pressure based design method in 
that the change in surface geometry is related to a change in pressure. The user starts with an initial 
"rough" cut of the geometry. The initial geometry is gridded and then analyzed using a CFD code. The 
surface pressure is then examined and a target pressure distribution is specified which is an 
improvement to the original pressure distribution. The difference between the desired pressure and the 
existing pressure is used to compute a change in geometry. The volume grid is automatically adjusted to 
accommodate the new geometry and a fresh analysis is performed. The results from this analysis are 
then compared to the target and a new geometry is computed and analyzed. The process continues until 
the user is satisfied with the analysis pressure. 

The nature of the relationship between geometry and pressure is different for subsonic (where change in 
pressure is related to change in curvature) and supersonic (where change in pressure is related to change 
in slope) flow but the basic idea is the same for both: that the resultant change in geometry depends on 
the desired change in pressure. This kind of design philosophy is ideally suited for nacelle design 
because nacelle flows are largely attached. Without special treatments, one would expect DISC to fail in 
regions of massive separation because the streamline curvature relationship would no longer be valid on 
the surface. 

In CDISC, one specifies desirable characteristics of the target pressure distribution. For example, one 
might specify pressure gradients in one region, shock position in another region, pressure levels in a 
third and then smooth selectively certain critical regions of the design space. These regions can be both 
chordwise and spanwise (circumferential) in extent. The target file is automatically generated and 
updated at each step based on all these constraints. One can also specify geometry constraints such as 
minimum thickness and leading edge radius. 

Additionally, CDISC has some global constraints for quantities l ik e twist distribution, lift coefficient 
and pitching moment that are especially useful for the design of wings. In the current version of CDISC 
an attempt is made to satisfy all constraints by iterating through them a set number of times in the order 
in which they are listed. 

Implementation 

DISC and CDISC have been coupled to many CFD codes. Airplane winglets and installed nacelles have 
been designed using DISC coupled to an Euler code4&5. DISC has been used with the CFL3D 
Navier-Stokes code for the design of isolated nacelles6. 

The coupling of CDISC with the multi-block codes, TLNS3D-MB (Thin-Layer Navier-Stokes 
3-Dimensional Multi-Block from Vatsa7,8&9) and OVERFLOW (a 3-D Navier Stokes code for overset 
Grids from BuninglO-13) is presented here. These two implementations of CDISC allow for the design 
of either isolated or installed nacelles with any combination of toe, inclination and roll angle for the 



nacelle. While this paper concentrates on nacelle design, it is important to note that the implementations 
are sufficiently general to allow for the redesign of any other part of the airplane. 

TLNS3D-MB 

CDISC is incorporated as a system of subroutines called by the main TLNS3D-MB driver. The CDISC 
routines are only called when the residual has dropped to a value chosen by the user. Once called, 

CDISC first extracts all relevant volume grids for the cowl surface from the grid array. These volume 
grids are then oriented to place the nacelle centerline on the x-axis. The volume grid is now on a sort of 
"virtual" workbench on which the design changes for the current design cycle can be easily performed. 

At each chosen streamwise design station, positioned circumferentially around the nacelle, the geometry 
is changed to offset the discrepancy between the analysis and the target pressures. The portions of the 
nacelle between target stations are changed in a linear fashion. Finally, circumferential smoothing is 
applied as desired. The volume grids are then perturbed to accommodate this new nacelle surface. 
Finally, the new grids are restored to their proper orientation and are placed back into the grid array for 
further analysis. This new analysis proceeds until a user specified drop in the residual between design 
iterations is achieved, after which a new design cycle is begun. 

The TLNS3D-MB solver, used herein, incorporates mesh sequencing and multigrid, but requires 
point-matched grids. The code allows a choice of turbulence models. The Spalart-Allmaras one-equation 
modell4 was used for the examples in this paper. 

One limitation in the above approach is the difficulty in constructing point-matched multi-block grids 
for complex configurations. For this reason, the overset grid approach is more attractive although there 
might be reservations among some in the aerospace community regarding its accuracy, particularly in 
regard to flux conservation. It should be noted, however, that most patched grid multi-block codes use 
non-conservative interpolation for block interfaces as well. 

OVERFLOW 

Except for the process of monitoring when to start and stop analysis and design, the entire CDISC 
system is largely run as a post-processor outside OVERFLOW. Control passes out of OVERFLOW 
whenever a user specified criterion (typically, the residual drop on the grid involving the component to 
be designed) is reached. Once out of OVERFLOW, the volume grid and restart files are used to compute 
the pressure on the fan cowl. The relevant grids are then oriented on to the "virtual" workbench 
described above and the CDISC geometry modifier is invoked. The grid is then automatically adjusted 
to accommodate the geometry change. The cowl grids are then restored to their proper orientation. 

Next, all relevant grids are re -projected and recomputed. For example, the pylon might have to be 
reprojected back onto the new nacelle surface. Next the pylon grid will have to be restretched and the 
nacelle fringe points re-projected back onto the pylon grid. Finally, the new grid system must be run 
through the mesh interpolation program, PEGSUS15, in order to update the interpolation stencils 
between the various grids, before OVERFLOW can be restarted for analysis of the new geometry. The 
analysis continues until the next residual drop is achieved. 

The OVERFLOW code also allows a choice of turbulence models. The Baldwin-Barth modell6 was 
used for the examples given in this paper. 


Results 



Isolated nacelle 

An example of designing an isolated nacelle using TLNS3D-MB is shown in Figures 1-4. The 
geometry is an isolated fan-cowl and core cowl. The grid for the fan-cowl and core-cowl comprised six 
blocks containing a total of one million grid points. The core cowl was run inviscid while the fan cowl 
was run viscous. The external surface of the fan cowl was redesigned. 

Figure 1 shows the convergence history for the design of this isolated nacelle. The spikes indicate the 
design iterations. The horizontal axis on the plot is a measure of the computational work that was 
performed. It is readily seen that in order to achieve the same 0(3.5) drop in residual, the design takes 
roughly 420/290 ~ 1.45 times the computational work of an analysis of a fixed geometry. Similar run 
comparisons were observed with the OVERFLOW runs. 

Figure 2 shows that the target coefficient of pressure has been closely matched after eight design cycles. 
The target stations were grid lines closest to the 45 degrees, 90 degrees, 135 degrees, 180 degrees, 225 
degrees, 270 degrees and 315 degrees cuts respectively. The run was set up so that no design changes 
were allowed on the upper part of the nacelle between 340 degrees and 20 degrees. One can see that the 
final Cp distribution is quite different from the initial distribution. 

This case is a good check-out run because the target pressure distribution was deliberately obtained from 
a geometry that was known a priori. While there is no guarantee that the inverse design will produce this 
identical geometry, Figures 3 and 4 show that the designed geometry is quite close to the "target" 
geometry. The dotted lines for the design shape are almost completely masked by the solid lines for the 
known "target" geometry. Note also that the resultant nacelle is symmetrical about the vertical plane. 
This symmetry condition was not forced during the design by using the CDISC option of aliasing 
corresponding left side stations to those on the right. Rather, the symmetry of the original geometry was 
maintained because of the uniform inflow condition. Note also that the original and designed nacelles 
are not symmetrical, about a horizontal plane passing through the nacelle centerline, because of the inlet 
droop. 

This example shows that a desired nacelle can be obtained from a radically different initial profile. The 
technique is directly applicable to any multi-block grid containing an installed nacelle. 

Installed Nacelle 

The nacelle designed above was installed on a transport airplane operating at a mid-cruise Mach number 
of 0.85. This flow-through nacelle was installed with a toe-in of 1.5 degrees and an upward tilt of 2 
degrees. The cruise angle of attack was roughly 1.5 degrees which meant that the installed nacelle faced 
an inflow stream at a much higher effective angle of attack than in the example above. Also the inflow 
was non-uniform because of the presence of the airplane wing and fuselage. 

This redesign was performed with OVERFLOW-CDISC. The fuselage-wing, pylon, fan-cowl, 
core-cowl, bifurcator, pylon shelf and far- field are represented by seven overset grids. The grid topology 
is not presented here because it is very similar to that described in reference 10 which also shows 
various views of the overset grid system. During each design cycle, the fan-cowl grid was the only grid 
extracted to the virtual workbench. The design was carried out over the entire circumference of the fan 
cowl. 


Figure 5 is a plot of the convergence history. The nine spikes indicate the design iterations. An 
erroneous tenth spike is also visible. This spike may be ignored; it occurred because of a malfunction in 



the computer system. Analysis of the final design change was stopped after 0(3.5) convergence was 
achieved. 


Figure 6 shows the Cp and geometry history for the nine design cycles at an inboard cut, the keel line, 
and an outboard cut. In this example, left-right symmetry is not maintained. The design is an attempt to 
achieve a favourable pressure gradient over the first 50 percent of the nacelle chord, while still allowing 
the inboard side of the nacelle to have overall lower Cp levels than the outboard side. This concession is 
not necessary, but is typical of the Mach suppression effect seen on the inboard side of an installed 
nacellel7. With CDISC one was able to achieve the characteristics of a desired pressure distribution 
without compromising the levels of the existing analysis pressure. 

Figure 7 is representative of the wing Cp, which was largely unaffected by the design changes to the 
nacelle. Although not shown, the wing drag did not change despite the dramatic changes in nacelle Cp. 

Conclusions 

Tools have been developed for the inverse design of a component of a complete airplane configuration 
while it is in the presence and under the influence of the rest of the airplane. Except for the initial grid 
generation, which will depend on the particular problem of interest, the design process is automated. 

The results showed the specific application of these tools to the design of turbofan nacelles. 
Additionally, it was shown that computational cost of the design is of the same order as that of a 
comparable analysis. 
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Fig. 1 Convergence history For inverse design of isloated raccllc using TLNS3D WB DISC. The peak 
residual (not visible in the piot) was 2.0S on thee oars estmeash. Two hundred cycles ofmesh 
sequencing were performed on two coarser levds beforethrer? level mulbgrid was started. 
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Fi* 2 Target. initial and design Cp For isolated nacelle. After eight design cycles. 
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